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Abstract. Entanglement is a non local property of quantum states 
which has no classical counterpart and plays a decisive role in quantum 
information theory. Several protocols, like the teleportation, are based 
on quantum entangled states. Moreover, any quantum algorithm which 
does not create entanglement can be efficiently simulated on a classical 
computer. The exact role of the entanglement is nevertheless not well un- 
derstood. Since an exact analysis of entanglement evolution induces an 
exponential slowdown, we consider approximative analysis based on the 
framework of abstract interpretation. In this paper, a concrete quantum 
semantics based on superoperators is associated with a simple quantum 
programming language. The representation of entanglement, i.e. the de- 
sign of the abstract domain is a key issue. A representation of entangle- 
ment as a partition of the memory is chosen. An abstract semantics is 
introduced, and the soundness of the approximation is proven. 



1 Introduction 

Quantum entanglement is a non local property of quantum mechanics. The en- 
tanglement reflects the ability of a quantum system composed of several sub- 
systems, to be in a state which cannot be decomposed into the states of the 
subsystems. Entanglement is one of the properties of quantum mechanics which 
caused Einstein and others to dislike the theory. In 1935, Einstein, Podolsky, 
and Rosen formulated the EPR paradox [7]. 

On the other hand, quantum mechanics has been highly successful in pro- 
ducing correct experimental predictions, and the strong correlations associated 
with the phenomenon of quantum entanglement have been observed indeed [5] . 

Entanglement leads to correlations between subsystems that can be exploited 
in information theory (e.g., teleportation scheme [3]). The entanglement plays 
also a decisive, but not yet well-understood, role in quantum computation, since 
any quantum algorithm can be efficiently simulated on a classical computer when 
the quantum memory is not entangled during all the computation. As a conse- 
quence, interesting quantum algorithms, like Shor's algorithm for factorisation 
|19| . exploit this phenomenon. 

In order to know what is the amount of entanglement of a quantum state, 
several measures of entanglement have been introduced (see for instance [13]). 
Recent works consist in characterising, in the framework of the one-way quantum 



computation |20| . the amount of entanglement necessary for a universal model 
of quantum computation. Notice that all these techniques consist in analysing 
the entanglement of a given state, starting with its mathematical description. 

In this paper, the entanglement evolution during the computation is analysed. 
The description of quantum evolutions is done via a simple quantum program- 
ming language. The development of such quantum programming languages is 
recent, see |17l8j for a survey on this topic. 

An exact analysis of entanglement evolution induces an exponential slow- 
down of the computation. Model checking techniques have been introduced [9] 
including entanglement. Exponential slowdown of such analysis is avoided by 
reducing the domain to stabiliser states (i.e. a subset of quantum states that can 
be efficiently simulated on a classical computer). As a consequence, any quantum 
program that cannot be efficiently simulated on a classical computer cannot be 
analysed. 

Prost and Zerrari |16| have recently introduced a logical entanglement anal- 
ysis for functional languages. This logical framework allows analysis of higher- 
order functions, but does not provide any static analysis for the quantum pro- 
grams without annotation. Moreover, only pure quantum states are considered. 

In this paper, we introduce a novel approach of entanglement analysis based 
on the framework of abstract interpretation 0. A concrete quantum seman- 
tics based on superoperators is associated with a simple quantum programming 
language. The representation of entanglement, i.e. the design of the abstract 
domain is a key issue. A representation of entanglement as a partition of the 
memory is chosen. An abstract semantics is introduced, and the soundness of 
the approximation is proved. 

2 Basic Notions and Entanglement 
2.1 Quantum Computing 

We briefly recall the basic definitions of quantum computing; please refer to 
Nielsen and Chuang [13] for a complete introduction to the subject. 

The state of a quantum system can be described by a density matrix, i.e. a 
self adjoinlQ positive-semidefinit(|l complex matrix of trac(H less than one. The 
set of density matrices of dimension n is £)„ C C" ^ " . 

The basic unit of information in quantum computation is a quantum bit or 
qubit. The state of a single qubit is described by a 2 x 2 density matrix p £ D2. 
The state of a register composed of n qubits is a 2" x 2" density matrix. If two 
registers A and B are in states S I?2" and ps e D2^, the composed system 
A, B is in state pa® Pb € I?2"+'"- 

The basic operations on quantum states are unitary operations and measure- 
ments. A unitary operation maps an n-qubit state to an ri-qubit state, and is 

^ M is self adjoint (or Hermitian) if and only if = M 

^ M is positive-semidefinite if all the eigenvalues of M are non-negative. 

^ The trace of M (tr(M)) is the sum of the diagonal elements of M 



given by a 2" X 2"-unitary matrijil- If a system in state p evolves according to 
a unitary transformation [/, the resulting density matrix is UpW. The parallel 
composition of two unitary transformations Ua, Ub is Ua ^ Ub- 

The following unitary transformations form an approximative universal fam- 
ily of unitary transformations, i.e. any unitary transformation can be approxi- 
mated by composing the unitary transformations of the family |13| . 
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A measurement is described by a family of projectors {P^, x e X} satisfying 

= Pi, PiPj = if i 7^ j, and X^xeA' Px = I- A computational basis measure- 
ment is {Pfe,0 < k < 2"}, where Pfc has entries everywhere except one 1 at 
row fc, column k. The parallel composition of two measurements {P^tX € X}, 
{P;, 2/ e Y} is {P, ® P;, {x, y)eXx Y}. 

According to a probabilistic interpretation, a measurement according to 
{Fx,x G X} of a state p produces the classical outcome x G X with proba- 
bility tr(P2,pP^) and transforms p into tr(P^pP ) ^xpPx- 

Density matrices is a useful formalism for representing probability distri- 
butions of quantum states, since the state p of a system which is in state pi 
(resp. P2) with probability pi (resp. P2) is p = pipi +P2P2- As a consequence, a 
measurement according to {P^, x £ X} transforms p into X^^gx ^xP^x- 

Notice that the sequential compositions of two measurements (or of a mea- 
surement and a unitary transformation) is no more a measurement nor a unitary 
transformation, but a superopcrator, i.e. a trace-decreasing completely positiv^ 
linear map. Any quantum evolution can be described by a superopcrator. 

The ability to initialise any qubit in a given state po, to apply any unitary 
transformation from a universal family, and to perform a computational mea- 
surement are enough for simulating any superopcrator. 

2.2 Entanglement 

Quantum entanglement is a non local property which has no classical counter- 
part. Intuitively, a quantum state of a system composed of several subsystems is 

* [/ is unitary if and only ii UW ^ UW = I. 
F IS. trace decreasing iff ti{F(p)) < tr(p) for any p in the domain of F. Notice 
that superoperators are sometimes defined as trace-perserving maps, however trace- 
decreasing is more suitable in a semantical context, see [TS] for details. 

® F is positive if F(p) is positive-semidefinite for any positive p in the domain of F. 
F is completely positive if Ifc ® _F is positive for any k, where Ife : C'''^'" — » C'"'^'" is 
the identity map. 



entangled if it cannot be decomposed into the state of its subsystems. A quantum 
state which is not entangled is called separable. 

More precisely, for a given finite set of qubits Q, let n = \Q\. For a given 
partition A, B of Q, and a given p £ D2'>, P is biseparable according to A, B (or 
{A, B)-separable for short) if and only if there exist K, pk > 0, p^ and such 
that 



p is entangled according to the partition A, B if and only if p is not {A, B)- 
separable. 

Notice that biscparability provides a very partial information about the en- 
tanglement of a quantum state, for instance for a 3-qubit state p, which is 
({!}, {2, 3})-separable, qubit 2 and qubit 3 may be entangled or not. 

One way to generalise the biseparability is to consider that a quantum state 
is TT-separable - where n = {Qj,j S J} is a partition of Q - if and only if there 
exist K, Pk > 0, and p^^ such that 



Notice that the structure of quantum entanglement presents some interesting 
and non trivial properties. For instance there exist some 3-qubit states p such 
that p is bi-separable for any bi-partition of the 3 qubits, but not fully separa- 
ble i.e., separable according to the partition {{1}, {2}, {3}}. As a consequence, 
for a given quantum state, there is not necessary a best representation of its 
entanglement. 

2.3 Standard and diagonal basis 

For a given state p £ T)'^ and a given qubit g e Q, if p is ({g}, (5\ {q})-scparable, 
then q is separated from the rest of the memory. Moreover, such a qubit may 
be a basis state in the standard basis (s) or the diagonal basis (d), meaning 
that the state of this qubit can be seen as a 'classical state' according to the 
corresponding basis. 

More formally, a qubit g of p is in the standard basis if there exists poiPi > 0, 
and po, Pi G g^ch that p = poP*'"^ ® Po+PiPl^^^^ ® Pi- Equivalcntly, q is 

in the standard basis if and only if P^'"^ pPl"^^^ = P^^^^^pP^'"^ = 0. A qubit q is 
in the diagonal basis in p if and only if q is in the standard basis in HqpHq. 

Notice that some states, like the maximally mixed 1-qubit state i(P'™^ -|- 
pfaise-j g^j.g both standard and diagonal basis, while others are neither in stan- 
dard nor diagonal basis like the 1-qubit state THP'^'"^HT . 

We introduce a function (3 : — > B'^ , where B'^ = Q ^ {s, d, T, _L}, such 
that /3(p) describes which qubits of p are in the standard or diagonal basis: 




Definition 1. For any finite Q, let P : V'^ — > B'^ such that for any p e 
and any q (z Q, 

if q is in both standard and diagonal basis in p 
if q is in the standard and not in the diagonal basis in p 
if q is in the diagonal and not in the standard basis in p 
otherwise 

3 A Quantum Programming Language 

Several quantum programming languages have been introduced recently. For a 
complete overview see [5] . We use an imperative quantum programming language 
introduced in [15| , the syntax is similar to the language introduced by Abramsky 
[1]. For the sake of simplicity and in order to focus on entanglement analysis, the 
memory is supposed to be fixed and finite. Moreover, the memory is supposed to 
be composed of qubits only, whereas hybrid memories composed of classical and 
quantum parts are often considered. However, contrary to the quantum circuit 
or quantum Turing machine frameworks, the absence of classical memory does 
not avoid the classical control of the quantum computation since classically- 
controlled conditional structures arc allowed (sec section [3TT1 ) 

Definition 2 (Syntax). For a given finite set of symbols q Cz Q, a program is 
a pair (C, Q) where C is a command defined as follows: 

C ::= skip 

I C*! ; 6*2 

I if q then Ci else C2 
I while q do C 

|H(g) 
\T{q) 

I CNot(g,q) 

Example 1. Quantum entanglement between two qubits 92 and 93 can be created 
for instance by applying H and CNot on an appropriate state. Such an entangled 
state can then be used to teleporte the state of a third qubit qi . The protocol of 
teleportation [2j can be described as (teleportation, {gi, 52, 93}), where 

teleportation : H{q2); 

CNot(g2,g3); 
CNot(gi,g2); 
H(gi); 

if qi then 

if q2 then skip else crx(<Z3) 
else 

if q2 then aziqs) else ayiq^) 
The semantics of this program is given in example [5) 




3.1 Concrete Semantics 



Several domains for quantum computation have been introduced |12I1I14| . 
Among them, the domain of superoperators over density matrices, introduced 
by Sehnger |18| turns out to be one of the most adapted to quantum semantics. 
Thus, we introduce a denotational semantics following the work of Sclinger. 

For a finite set of variables Q = {qo, . . . , qn}, let T)'^ = D2\q\ . Q is a set of 
qubits, the state of Q is a density operator in V'^ . 

Definition 3 (Lowner partial order). For matrices M and N m C"^", M C 
N if N ~ M is positive-semide finite. 

In [TB], Selinger proved that the poset iV^ ,'^) is a complete partial order 
with as its least element. Moreover the poset of superoperators over IS a 
complete partial order as well, with as least element and where the partial order 
a is defined as F C' G Vfc > 0, Vp e V^^^ , (Ife ® F){p) C (1^ G)(p), 

where Ifc : Vk is the identity map. Notice that these complete partial 

orders are not lattices (see [HI-) 

We are now ready to introduce the concrete denotational semantics which 
associates with any program (G, Q) , a superoperator |G] : V'^ . 

Definition 4 (Denotational semantics). 

[skip] = I 

ICi;G2] = IG2]o[Gil 
[U(g)l - Xp.U.pUl 
[CNot(qi,92)l = \p.CNotq,,q,pCNot\^^^^ 

lif q then Gi else G2 ] = \p. ([GiKPrVP'r) + IC2l(P;='=VP^^'^)) 

[while g do G 1 = Ifp (\.f.\p. (/ o [GMP^VPr) + P'fVP'f 0) 

where P'™^ " (o o) ^'"''^ " (o l)' " \p.MpM\ and Mq means 
that M is applied on qubit q. We refer the reader to an extended version of this 
paper for the technical explanations on continuity and convergence. 

In the absence of classical memory, the classical control is encoded into the 
conditional structure if q then Gi else G2 such that the qubit q is first measured 
according to the computational basis. If the first projector is applied, then the 
classical outcome is interpreted as true and the command Gi is applied. Other- 
wise, the second projector is applied, and the command G2 is performed. The 
classical control appears in the loop while g do G as well. 



As a consequence of the classical control, non unitary transformations can 
be implemented: 



|if q then q else a. 



X 



(q) j : P{«> ^ P{«> = Ap.P' 



itrue 



[while q do H(g) ] : V^i^ -> V^i^ = Xp.P 



I false 



Notice that the matrices P*™^ and P ^ used in dcfinition[3]for describing the 
computational measurement {pt''"^ pf^'^^j. can also be used as density matrices 
for describing a quantum state as above. 

Moreover, notice that all the ingredients for approximating any superopera- 
tors can be encoded into the language: the ability to initialise any qubit in a given 
state (for instance p'™^ or P*^^'^^); an approximative universal family of unitary 
transformation {H, T, CNot, a^, Uy, cr^}; and the computational measurement of 
a qubit q with if q then skip else skip . 

Example 2. The program (teleportation, {qi, (72, 93}) described in example [T] re- 
alises the teleportation from qi to 53, when the qubits 52 and 53 are both ini- 
tialised in state P*'^"^: for any p E T>2, 



4 Entanglement Analysis 

What is the role of the entanglement in quantum information theory? How 
does the entanglement evolve during a quantum computation? Wc consider the 
problem of analysing the entanglement evolution on a classical computer, since 
no large scale quantum computer is available at the moment. Entanglement 
analysis using a quantum computer is left to further investigation^ 

In the absence of quantum computer, an obvious solution consists in sim- 
ulating the quantum computation on a classical computer. Unfortunately, the 
classical memory required for the simulation is exponentially large in the size of 
the quantum memory of the program simulated. Moreover, the problem SEP of 
deciding whether a given quantum state p is biseparablc or not is NP HarcJlUn]- 
Furthermore, the input of the problem SEP is a density matrix, which size is 
exponential in the number of qubits. As a consequence, the solution of a classical 
simulation is not suitable for an efficient entanglement analysis. 

^ Notice that this is not clear that the use of a quantum computer avoids the use of the 
classical computer since there is no way to measure the entanglement of a quantum 
state without transforming the state. 

® For pure quantum states (i.e. tr(p^) = tr(p)), a linear algorithm have been 
introduced [11] to solve the sub-problem of finding biseparability of the form 
({go, ■ ■ ■ ,<lk}, {qk+i , ■ ■ ■ ,<ln}) - thus sensitive to the ordering of the qubits in the 
register. Notice that this algorithm is linear in the size of the input which is a den- 
sity matrix, thus the algorithm is exponential in the number of qubits. 




To tackle this problem, a solution consists in reducing the size of the quan- 
tum state space by considering a subspace of possible states, such that there 
exist algorithms to decide whether a state of the subspace is entangled or not in 
a polynomial time in the number of qubits. This solution has been developed in 
[9], by considering stabiliser states only. However, this solution, which may be 
suitable for some quantum protocols, is questionable for analysing quantum algo- 
rithms since all the quantum programs on which such an entanglement analysis 
can be driven are also efficiently simulablc on a classical computer. 

In this paper, we introduce a novel approach which consists in approximating 
the entanglement evolution of the quantum memory. This solution is based on 
the framework of abstract interpretation introduced by Cousot and Cousot [5]. 
Since a classical domain for driving a sound and complete analysis of entangle- 
ment is exponentially large in the number n of qubits, we consider an abstract 
domain of size n and we introduce an abstract semantics which leads to a sound 
approximation of the entanglement evolution during the computation. 

4.1 Abstract semantics 

The entanglement of a quantum state can be represented as a partition of the 
qubits of the state (see section [52]) , thus a natural abstract domain is a domain 
composed of partitions. Moreover, for a given state p, one can add a flag for each 
qubit q, indicating whether the state of this qubit is in the standard basis s or 
in the diagonal basis d (see section [2?3| . 

Definition 5 (Abstract Domain). For a finite set of variables Q, let A"^ = 
qQ X Y\Q be an abstract domain, where = Q ^ {s,d,T,±} and fl'^ is the 
set of partitions of Q: 

n<5 = {tt C p(Q) \ {0} I [j X^Q and (VX, F e tt, X C]Y = ^ or X = Y)} 

The abstract domain A is ordered as follows. First, let ({s,d,T,_L},<) be a 
poset, where < is defined as: _L < s < T and _L < d < T. {B^ , <) is a poset, 
where < is defined pointwisc. Moreover, for any 7ri,7r2 G . let tti < 112 if i^i 
rafines 7r2, i.e. for every block X € tti there exists a block Y <^ 112 such that 
X <ZY. Finally, for any (61, tt), (62, 772) £ A'^ , (&i,7r) < (62,7r2) if 61 < 62 and 

TTi < 7r2. 

Proposition 1. For any finite set Q, {A'^ , <) is a complete partial order, with 
_L = (Aq.-L, {{q}, q £ Q}) as least element. 

Proof. Every chain has a supremum since Q is finite. □ 

Basic operations of meet and join are defined on A'^ . It turns out that con- 
trary to V^, {A'^,\/,A,±,{Xq.T,{Q})) is a lattice. 

A removal operation on partitions is introduced as follows: for a given parti- 
tion TT = {Qi, « G /}, let TT \ g = {Qi \ {q}, i € /} U {{(?}}. Moreover, for any pair 
of qubits gi, 52 € Q, let [gi, 92] = {{q \ q & Q \ {<?!, 92}}, {gi, 92}}- 



Finally, for any b G B'^ , any qo,q & Q, any k £ {s, d, T, ±}, let 



k if g = go 
bq otherwise 



We are now ready to define the abstract semantics of the language: 

Definition 6 (Denotational abstract semantics). For any program (C, Q), 
let ICf ■.A'^^A'^ be defined as follows: For any (6, ir) (E A^ , 



IskiplH^,^) 



(fe,7r) 

lC2folC,f{b,7r) 



laiq)fib,TT) = {b,n) 



lH{q)f{b,7r) = 



[T(<z)f(6,^) 



ICNot(gi,g2)lH^^) 



W^'^,n) ifbg = s 
tt) ifbq^d 
b, tt) otherwise 

b-i"^^ ,t:) ifbq^d 

&, tt) otherwise 



b, tt) 



|if q then Ci else C2 f{b,Ti) 
[while g do C 1^6, tt) 



i/ = s or = d 

tt) i/ = ± and > ± 

tt) i/ > ± and bq.^ = ± 
^giH^s,92^d^ tt) = ± and bg^ = ± 

fe'?!'''^'^^, TT V [qi, (72]) otherwise 
(ICif (&'^^% ^ \ V [C2l^6«^^ TT \ g)) 
Ifp (A/.A^. (/ o ICfib'i-^, 7T\q)V TT \ g))) 



where F^ = A(6, 7^).(6«'-*^ tt \ q). 

Intuitively, quantum operations act on entanglement as follows: 

— A 1-qubit measurement makes the measured qubit separable from the rest 
of the memory. Moreover, the state of the measured qubit is in the standard 
basis. 

— A 1-qubit unitary transformation does not modify entanglement. Any Pauli 
operator a G {(Jx, cTy, az} preserves the standard and the diagonal basis of 
the qubits. Hadamard H transforms a state of the standard basis into a 
state of the diagonal basis and vice-versa. Finally the phase T preserves the 
standard basis but not the diagonal basis. 



— The 2-qubit unitary transformation CNot, applied on qi and q2 may create 
entanglement between the qubits or not. It turns out that if qi is in the stan- 
dard basis, or q2 is in the diagonal basis, then no entanglement is created and 
the basis of qi and q2 are preserved. Otherwise, since a sound approximation 
is desired, CNot is abstracted into an operation which creates entanglement. 

Remark 1. Notice that the space needed to store a partition of n elements is 
0(n). Moreover, meet, join and removal and can be done in cither constant or 
linear time. 

Example 3. The abstract semantics of the tclcportation (see example [1]) is 
|teleportation|h : y^i-Ji'-J^,?,,} ^ _4{'?i,92,93} = A(6, 7r).(6«i'«^^^'''-^""^, _L). Thus, 
for any 3-qubit state, the state of the memory after the teleportation is fully 
separable. 

Assume that a fourth qubit 94 is entangled with qi before the telepor- 
tation, whereas q2 and (J3 are in the state P'™^. So that, the state of the 
memory before the teleportation is [gi, g4]-separable. The abstract semantics 
of (teleportation, {gi, 92, 93, qi]) is such that 

Iteleportationf (6, [91,94]) - (fe^i'^^-^^^-T^ 

Thus the abstract semantics predicts that (73 is entangled with 94 at the end 
of the teleportation, even if 93 never interacts with 94. 

Example 4. Consider the program (trap, {qi, (72}), where 
trap = CNot(9i,92);CNot(gi,92) 
Since CNot is self- inverse, |trap| : X'^'Ji'*^} _^ pin, 92} = Xp p Yoi instance, 

|trap](i(Pt™^ ^ pfalse-) ^ ptrue-) ^ l(ptrue _^ pfalse-j ^ ptrue^ 

However, if bq-^ ~ d and bq^ = s then 

[traplH6,{{<zi},{92}}) = (fo^^^^'^^^^, {{gi, 92}}) 

Thus, according to the abstract semantics, at the end of the computation, qi 
and q2 are entangled. 



4.2 Soundness 

Example [3] points out that the abstract semantics is an approximation, so it may 
differ from the entanglement evolution of the concrete semantics. However, in 
this section, we prove the soundness of the abstract interpretation (theorem [T]) . 

First, we define a function /3 : V'^ such that /3(p) describes which 

qubits of p are in the standard or diagonal basis: 

Definition 7. For any finite Q, let (3 : B'^ such that for any p G 

and any q <E Q, 

{s if P'^'^^pPl"''" = Pl''"'pP'g""' = 
d if (p^™= + pf="==)p(p^™= - Pi^'n = (p'r - p'tip{p'r + K^'n = 
T otherwise 



A natural soundness relation is then: 

Definition 8 (Soundness relation). For any finite set Q, let a G p(I?'^,^'^) 
he the soundness relation: 

a ~ {(p, (6,7r)) I p is n-separable and P{p) < 6} 

The approximation relation is nothing but the partial order <: (6, tt) is a more 
precise approximation than (6',7r') if (6, tt) < (6',7r'). Notice that the abstract 
soundness assumption is satisfied: if p is 7r-separable and tt < tt' then p is tt'- 
separable. So, (p, a) g cr and {p,a) < {p',a') imply {p',a') E a. 

However, the best approximation is not ensured. Indeed, there exist some 
3-qubit states |6l4j which are separable according to any of the 3 bipartitions of 
their qubits {a, b, c} but which are not {{a}, {6}, {c}}-separable. Thus, the best 
approximation does not exist. 

However, the soundness relation a satisfies the following lemma: 

Lemma 1. For any finite set Q, any pi,P2 G 'D'^ , and any ai,a2 G A'^ , 

(pi, ai), (/02, 02) e (T =^ {pi + P2, TTl V 7r2) e CT 

Moreover, the abstract semantics is monotonic according to the approxima- 
tion relation: 

Lemma 2. For any command C, fCJ^ is <-monotonic: for any tti,tt2 G A'^ , 

TTl <7r2 =^ [ClH^i) < [C7f(7r2) 
Proof. The proof is by induction on C. 

Theorem 1 (Soundness). For any program {C,Q), any p e V'^ , and any 

aeAQ, 

{p,a)ea =^ mip),lCf{a))ea 

Proof. The proof is by induction on C. 

In other words, if p is 7r-separable and (3{p) < b, then |C](/9) is 7r'-separable 
and mCKp)) < b', where (6',7r') = lCfib,7r). 

5 Conclusion and Perspectives 

In this paper, we have introduced the first quantum entanglement analysis based 
on abstract interpretation. Since a classical domain for driving a sound and 
complete analysis of entanglement is exponentially large in the number of qubits, 
an abstract domain based on partitions has been introduced. Moreover, since 
the concrete domain of superoperators is not a lattice, no Galois connection 
can be established between concrete and abstract domains. However, despite 



the absence of best abstraction, the soundness of the entanglement analysis has 
been proved. 

The abstract domain is not only composed of partitions of the memory, but 
also of descriptions of the qubits which arc in a basis state according to the stan- 
dard or diagonal basis. Thanks to this additional information, the entanglement 
analysis is more subtle than an analysis of interactions: the CNot transforma- 
tion is not an entangling operation if the first qubit is in the standard basis or 
if the second qubit is in the diagonal basis. 

A perspective, in order to reach a more precise entanglement analysis, is to 
introduce a more concrete abstract domain, adding for instance a third basis, 
since it is known that there are three mutually unbiased basis for each qubit. 

A simple quantum imperative language is considered in this paper. This 
language is expressive enough to encode any quantum evolution. However, a 
perspective is to develop such abstract interpretation in a more general setting 
allowing high-order functions, representation of classical variables, or unbounded 
quantum memory. The objective is also to provide a practical tool for analysing 
entanglement evolution of more sophisticated programs, like Shor's algorithm 
for factorisation [TO] . 

Another perspective is to consider that a quantum computer is available for 
driving the entanglement analysis. Notice that such an analysis of entanglement 
evolution is not trivial, even if a quantum computer is available, since a to- 
mography |21j is required to know the entanglement of the quantum memory 
stat^. 
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